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Abstract 

We examine in the context of supersymmetric models whether the usual 
cosmological upper bound on the axion decay constant can be relaxed by 
assuming a period of stronger QCD in the early universe. By evaluating the 
axion potential in the early universe and also taking into account the dilaton 
potential energy, it is argued that a stronger QCD is not useful for raising up 
the bound. 
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I. INTRODUCTION 



One of the most attractive solutions to the strong CP problem is to introduce a spon- 
taneously broken Peccei-Quinn symmetry This solution originally predicted a pseudo- 
Goldstone boson, the so-called Peccei-Quinn- Weinberg- Wilczek axion [[|, which is ruled out 
phenomenologically. The subsequently devised invisible axion 0] is also constrained to a 
narrow band on the axion decay constant f a . Some astrophysical arguments, e.g. those as- 
sociated with the axion emission from helium burning red giants or the supernova SN1987A 
H|J, imply that f a is far above the weak scale. A natural possibility is then that f a is 
around the grand unification scale or the Planck scale. As is well known, however, the co- 
herent axion oscillation in the early universe gives rise to relic axions whose mass density 
Q a (in the unit of the critical energy density) at present is given by || 

fia=(W/a) 2 (/a/4xl0 12 GeV) 1 - 18 , (1) 

where 5a is the axion misalignment at the epoch of the QCD phase transition. Clearly it 
is most natural to assume that the axion misalignment is of order f a . Requring the axion 
mass density not significantly exceed the critical density, this leads to the upper bound 
fa < 4 x 10 12 GeV which is far below the grand unification scale or the Planck scale. 

One way to raise up the cosmological upper bound on f a is to have an entropy production 
after the QCD phase transition @. However it has been argued that, due to constraints 
from the nucleosynthesis, one cannot reach above few times 10 15 GeV by this late time 
entropy production ||. Recently another interesting way to raise up the axion bound was 
suggested by Dvali |§. The idea is to suppress the axion misalignment 5a by assuming 
a period of stronger QCD in the early universe during which the QCD coupling constant 
takes a value larger than the present one. Clearly this would be possible only when the 
QCD coupling constant is determined by a dynamical degree of freedom in the theory, e.g. 
a dilaton field whose vacuum expectation value (VEV) gives 1/g 2 = (<p). If the dilaton 
takes a value smaller than the present one, the QCD at that time would be stronger, leading 
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to the corresponding scale Aq C d far bigger than the present value of order 0.2 GeV. Then 
there is a possibility that a large effective axion mass is induced in the early universe, driving 
the axion field toward the minimum of the effective potential. If the minimum of the early 
axion potential coincides with the minimum of the present axion potential, 5a would be 
dynamically suppressed, allowing the upper bound of /„ larger than 4 x 1 12 GeV. Note that 
Eq. (f) indicates that the upper bound of f a can be around 10 16 ~ I0 18 GeV if Sa/f a is 
relaxed down to ICT 2 ~ ICT 3 in the early universe. 

In order for the axion misalignment 5a to be relaxed by the above mechanism, one needs 

(i) m a >H, (ii) \(a/f a )-6 cS \<10- 3 ~10~ 2 , (2) 

where m a and (a) denote the axion mass and the vacuum expectation value (VEV) deter- 
mined by the axion potential in the early universe, and 9 e g is the low energy QCD vacuum 
angle at the present. Here the first condition is required for the onset of the axion field 
oscillation toward the minimum, and the second is for the location of this minimum to co- 
incide with the present one up to a small misalignment angle of order ICT 3 ~ fCT 2 . Note 
that the present axion VEV is given by the low energy vacuum angle 6 c k which includes the 
contributions from a variety of CP-odd parameters in the underlying theory. For instance, 
if the underlying theory is the minimal supersymmetric standard model, 

#cfr = 9 QC D + arg(det(A u A d )) + 3arg(mi /2 ) - 3arg(/i J B), (3) 

where Qqcd denotes the bare QCD vacuum angle, \ u and A^ are the Yukawa coupling 
matrices, and m^, /i and B denote the gluino mass, the Higgs fi parameter, and the soft 
B parameter, respectively. We stress that the second condition in Eq. (2) for the early 
minimum to coincide with the present one should be satisfied without fine tuning of some 
field values which would be chaotic in the early universe. Such fine tuning is certainly 
nothing better than the simpler fine tuning of the chaotic axion field a/ f a to 9 cS at the 
epoch of QCD phase transition. 

In this paper, we wish to examine whether the axion misalgnment can be relaxed by a 
stronger QCD in the early universe in supersymmetric models. To start with, we assume 
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that the dilaton dynamics in the model allows a period in the early universe during which 
QCD becomes stronger. The reason for limiting the analysis to supersymmetric models is 
clear: the most natural playground of the dilaton field is supersymmetric model. 

The organization of this paper is as follows. In Section 2, we provide a discussion of the 
axion potential in the early universe by distinguishing the two cases, one in which Aq C d is 
smaller than the size of soft supersymmetry breaking parameter m so ft and the other opposite 
case that the QCD becomes so strong that Aq C d 3> m soft . We discuss the second case in 
somewhat detail since the axion potential in this case has not been considered in the previous 
literature. As we will see, the axion potential in the second case dramatically differs from the 
well known axion potential in the first case. Based on this analysis, in Section 3 we examine 
whether the two conditions given in Eq. (2) can be simultaneously met. Taking into account 
the dilaton potential energy, we will argue that it is very unlikely to have a cosmological 
scenario in which both of these two conditions are satisfied. We thus conclude in Section 4 
that a stronger QCD in the early universe is not useful for raising the cosmological upper 
bound of the axion decay constant /„. 

II. AXION POTENTIAL IN THE PERIOD OF STRONGER QCD 

In this section, we present a somewhat detailed discussion of the axion potential in 
the early universe. In the early universe, field variables would take chaotic values if their 
effective masses are smaller than the expansion rate H. In order to satisfy Eq. (2) without 
fine tuning, one needs to keep the axion VEV not shifted by such chaotic fields. If not, 
the axion VEV in the early universe would become chaotic and thus not coincide with the 
present axion VEV in general. We will thus assume that field variables which apparently 
affect the axion VEV are not chaotic, but quickly settle into the minimum of their effective 
potential in the early universe. Of course, one still needs to confirm that these fields settled 
at the minimum leads to an early axion VEV almost the same as the present one as given 
in Eq. (2). 



4 



To compute the axion potential in the early universe, it is convenient to distinguish the 
following two cases: 

(?) Aqcd < m sott , (ii) Aqcd > m soft . 

Let us first consider the case that Aq C d < m aoit- In this case, if there is no light quark 
whose current mass is less than Aqcd, the axion potential is given by 

V* - A qcd cos(a// a - in ) , (4) 

where 9- in denotes the QCD vacuum angle in the early universe. The presence of light quarks 
suppresses the axion potential as 

V a ~ m q A 3 QCD cos(a// a - 9 in ), (5) 

where m q (< Aqcd) is the current mass of the lightest quark. 

Of course, it depends upon the Higgs VEVs in the early universe whether there exists a 
light quark with m q < Aqcd- In fact, the Higgs VEVs in the early universe can dramatically 
differ from the present ones. To be explicit, let us consider the minimal supersymmetric 
standard model. In this model, the Higgs potential for the neutral components takes the 
form 

m\\H u \ 2 + ml\H d \ 2 + (B»H v H d + h.c.) + l -(g\ + gl){\H u \ 2 - \H d \ 2 ) 2 , (6) 

where H u and denote the two Higgs doublets in the model. In generic supersymmetric 
models, soft parameters are determined by the VEVs of fields triggering spontaneous SUSY 
breaking. As a result, soft parameters in the early universe may take values which differ 
from the present ones. In particular, both the sign and magnitude of m\ and m\ in the early 
unverse can differ from the present ones due to an additional SUSY breaking by nonzero 
energy density. This change of m\ and m\ does not alter arg(H u H d ) and thus not lead to 
a change of the early axion VEV. However it can dramatically change \H U \ and \H d \. For 
instance, if both m\ and m 2 are positive and large enough, the Higgs VEVs do vanish. Then 
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all quarks become massless and the axion potential vanishes up to a tiny contribution from 
small instantons. In the opposite case that both m\ and m\ are negative, the Higgs VEVs 
slide down toward the flat direction \H U \ = \Hd\ ^> wi S oft- In this case, we would not have 
any light quark and the axion potential is given by (4). Finally for m\ and m\ leading to 
the Higgs VEVs similar to the present ones, the axion potential is given by (5) with the up 
quark mass m q ~ 10 MeV. 

In the above, we have ignored a possible finite temperature effect on the axion poten- 
tial. At high temperature T > Aqcd, the axion potentials Eq. (4) and Eq.(5) are not valid 
anymore. However, compared to the low temperature result, the high temperature axion po- 
tential is suppressed by more powers of the light quark masses and also by exp(— Sir 2 /g 2 (T)) 
T0|| . Since what we wish to obtain is a bigger axion potential for a given value of Aqcd, we 



will not consider such a high temperature case anymore. 

In the above discussion, the Higgs fields were assumed to stay at the minimum of the 
potential. It has been argued that in generic supergravity models a nonzero energy density 
Pe ~ H 2 Mp l in the early universe affects soft parameters as 5m so f t = cH where H is 
the Hubble expansion rate and c is a dimensionless constant [IT]. The constant c may be 



small due to small couplings or to loop suppression. If c is small enough, it is possible 
that m so f t <C H . Such a large H may result in chaotic Higgs values in the early universe. 
However, in supersymmetric models with two Higgs doublets H u and H^, chaotic Higgs 
fields lead to a chaotic axion VEV which differs from the present axion VEV in general. 
Furthermore, for H > m so f t , we have m a <C H (Note that we have been discussing the case 
Aqcd < wisoft-) an d thus the first condition of Eq. (2) can not be satisfied also. We thus 
assume that H < m so f t and thus the Higgs fields quickly settle into their VEVs. 

Let us now consider the opposite limit that the QCD becomes so strong that Aq C d ^ 
m so f t . Then the axion potential is induced by nonperturbative supersymmetric QCD effects. 
Recent progress in understanding the dynamics of supersymmetric gauge theories ]TJ|] allows 
us to estimate the axion potential even in this case. 

As is well known, the axion potential is severely affected by gauge invariant condensates 



which break chiral symmetries in the theory. These condensates can be used to tie together 
fermion zero modes of the instanton amplitudes generating the axion potential. Phase 
degrees of freedom of such condensates then mix with the axion field as the rf mixes with 
the axion in normal QCD J3j. We thus study first what kind of chiral symmetry breaking 
condensates are formed in the supersymmetric QCD with the superpotential 

W = \ u H u Qu c + \ d H d Qd c } (7) 

where Q, u c , and d c denote SU(2)l doublet quarks and singlet antiquarks, respectively. 
Here, the generation indices and gauge group indices are omitted for simplicity. 

It has been argued that, for a vanishing superpotential, the quantum moduli space of 
degenerate vacua for supersymmetric SU(3) C gauge theory with Nf = 6 quark flavors is the 



same as the classical one [13]. Such a vacuum degeneracy is lifted in fact by the Yukawa 



terms in the superpotential and also by soft breaking terms ||14j| . As was observed recently, 
the infrared behavior of gauge invariant operators in supersymmetric SU(N C ) gauge theory 
with Nf quark flavors can be studied by the dual SU (Nf — N c ) theory with Nf dual quarks 
]T3]] . The superpotential of the dual model contains 

W D = K QC d{KH u T u + \ d H d T d ) + T u Q D u c D + T d Q D d c D , (8) 

where T u = Qu c /Aq C d and T d = Qd c / Kq C d denote the meson superfields with dimension 
one at the ultraviolet, while Q D , u c D , and d c D are the dual quark and anti-quark superfields. 
Once the supersymmetry of the original model is broken by soft parameters much less than 



^qcdi the dual model will contain soft terms |14 

C^^AWn + Y.mWhW (9) 

where <pi denotes generic scalar fields in the dual model and the soft parameters A and mi 
are again much less than Aqcd- 

A nice feature of the dual theory is that it becomes weaker as the original QCD becomes 
stronger, allowing a classical approximation to be valid. The classical effective potential 



of the dual model can be readily computed from the superpotential (8) and the soft terms 
(9). It is then easy to see that, for Aq C d large enough compared to m soft , the Higgs fields 
H Ut d and the linear combinations \ U T U and XdTd of mesons have positive mass squared of 
order Aq CD , and thus have vanishing VEV's. For the original squarks having positive soft 
mass squared, it is expected that both the mesons and the dual squarks also have positive 
soft mass squared. Note that the mesons are the bound states of the original squarks and 
anti-squarks, while the dual squarks can be obtained by dissociating the scalar baryons 
(containing N c original squarks) into (Nf — N c ) pieces. With positive soft mass squared, the 
VEVs of the dual squarks and the entire mesons vanish also. 
The above results on the VEVs can be summarized as 

(qq c ) oc (T) = 0, 
(qq c ) oc (F T ) = 0, 
(AA) oc (T)JV(^/-^) = 0> (10 ) 

where T denotes the scalar components of the meson superfields T Utd with the auxiliary 
components F T . Thus all of the squark condensates (qq c ), the quark condensates (qq c ), and 
the gluino condensate (AA) do vanish. Note that the vanishing of the gluino condensate 
is essentially due to the vanishing Higgs VEV leading to the massless quarks. Although 
the above chiral invariant vacuum configurations were derived based on the classical poten- 
tial at zero temperature, it is rather easy to see neither quantum corrections (due to the 
weak dual gauge interactions) nor finite temperature effects does shift these chiral invariant 
configurations. 

So far, we have argued that the mesons and the dual squarks have vanishing VEVs for 
Aqcd m soit- Again if the expansion rate H is large enough so that H ^> m soft , the mesons 
and dual squarks would have chaotic values although their effective potential is minimal at 
zero values. Note that the dual squarks and the mesons (except for \ U T U and XdTd) have 
masses of order m so f t . Chaotic values of the mesons and the dual squarks correspond to 

8 



chaotic condensates of the quarks, squarks, and gluino. Clearly the axion VEV determined 
by such chaotic condensates will be chaotic also, and thus does not coincide in general with 
the present axion VEV. Thus to avoid a chaotic axion VEV, we assume that m so ft > H and 
the mesons and the dual squarks quickly settle into the minimum of their effective potential. 

In the case that there is no condensate, the axion potential is given by the QCD instanton 
amplitudes whose fermion zero modes are tied together by the interactions which break the 
chiral symmetries explicitly fl5fl . To proceed, let us consider the minimal supersymmetric 
standard model having the superpotential 

W = \ u H u Qu c + X d H d Qd c + fiH u H d} (II) 

and the soft terms 

Aoft = ^m 1/2 \\ + A(\ u H u Qu c + \ d H d Qd c ) + B^H u H d + i £>i l^l" + h.c, (12) 

where fa denotes generic scalar fields in the model. In order to compute the instanton- 
induced axion potential, we first note that the model is invariant under 

Gmssm = SU(3) Q x SU(3) u c x SU(3) dc x U(1) A x U(l) x x U(1) R , (13) 

where the fields and parameters transform as Q = (3, 1, 1), u c = (1, 3, 1), d c = (1, 1, 3), X u — 
(3,3, 1), and = (3, 1,3) under SU(3)q x SU(3) u c x SU(3) d c, and the quantum numbers 
of U(1)a x U(l)x x U(1)r are given in Table 1. As is well known, the soft supersymmetry 
breaking can be described within superspace formalism by introducing spurion superfields: 

v = {l+m 4 W}, 

Y = (1 + lQir 2 m 1/2 6 2 )T, 

Z = {Z u4 = (1 + A6 2 )\ u4l Z„ = (1 + Be 2 )^}, 

where the auxiliary components of the spurion superfields represent soft breaking, while 
the scalar components denote supersymmetric couplings: the complex gauge coupling r = 
^jr + i^QCD, the Yukawa couplings \ U:d , and the /i parameter. Note that the factor 16n 2 in 
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Y is introduced for a proper normalization of the gluino mass m^. Obviously r/s are real 
superfields while Y and Z's are chiral. 

After integrating out the gauge and matter superfields, the axion effective Lagrangian 
can be read off from the effective lagrangian of the spurions by identifying the dimensionless 
axion superfield 

A = {s + ia + a6 + F A 6 2 )/f a 
as the fluctuation of Y, i.e. by the identification 

Y^Y + A. (14) 

The effective Lagrangian of spurions can be written as 

J d 2 6d 2 6 K cS (Y, Y*, Z, Z*, 77) + J d 2 6 W cS (Y, Z) + h.c, (15) 

where the effective Kahler potential K e R is a real function of spurion superfields (and also 
of their supercovariant derivatives) while the effective superpotential W e g is a holomorphic 
function of chiral spurions. 

Since the matter and gauge fields are integrated over a unique ground state preserving 
chiral symmetries, the above effective Lagrangian does not have any branch cut associated 
with the multiplicity of ground states. Note that in the case with nonvanishing gaugino con- 
densate (AA) ~ e ~ ie QCD/N ^ corresponding effective Lagrangian contains terms depending 
upon e~ Y l N ~ e ~ te QCD/N anc j thug has N branches associated with the N degenerate vacua 



which are related to each other by the 2ir shift of Qqcd [ 16fl - In our case with a unique chiral 
invariant vaccum, the effective Lagrangian is a single valued function of e~ tdQCD , and thus 
is manifestly periodic under the 2ti shift of Qqcd- 
Instantons would induce a term in W e g as 

e~ nY uj{Z), (16) 

where n is a positive integer corresponding to the instanton winding number. It is, however, 
easy to see that the selection rules of Gmssm does not allow any holomorphic uj which is 
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finite at // — > 0. This implies that the axion potential does not appear through W e g, but 
through K e $. Again the selection rules of Gmssm implies that at the leading order the 
instanton-induced Kahler potential takes the form (for n — 1): 

K cS cx e- Y det(Z u Z d )Z; 3 F( V ) + h.c., (17) 

where F is an arbitrary function of the real superfields r\. 

To obtain the axion potential from the superspace integration of K e g, we need at least 
either a single insertion of D v or the simultaneous insertions of Fy,z and F Y z where D and F 
denote the auxiliary components of spurions. In other words, supersymmetry (together with 
the selection rules of Gmssm) implies that the axion potential is suppressed at least by the 
two powers of m so f t = {m,, m^, A, B}. Note that here we consider instantons whose scale, 
i.e. the inverse of the instanton size p, is in the range between m so f t (or p) and the messenger 
scale M m above which soft breaking is not operative anymore. In the popular hidden sector 



models in which supersymmetry breaking is transmitted by supergravity interactions |17 



the messenger scale M m = Mp\, while it can be much lower in visible sector models ||TB 
Obviously instantons with p < M" 1 (if exist) do not contribute to the axion potential due 
to the restored SUSY. As we will see, instantons which give dominant contributions have 
size p ~ ^-qcd an d thus belong to the above category for p ~ m so f t Aqcd- At any rate, 
from the Kahler potential of Eq. (17), we readily find the SUSY suppression factor 

K oft ] 2 = {ml AB\ I6n 2 m 1/2 B*}, (18) 

where the factor 16tt 2 in front of m^j indicates that instanton graphs using m^j to tie 
together gluino zero modes contain one less loop compared to those using other soft param- 
eters. 

With the above observation, one can write the axion potential in the case Aqcd 3> ^soft 

as 

V a = e ia/fa p* 3 det(\ u \ d ) tt + h.c. (19) 
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where Q is suppressed by [m so f t ] 2 and also by some powers of the loop factor 1/167T 2 . One 
may estimate Q using an explicit instanton graph. For instance, the dimensional analysis 
rule of Ref. [[HJ applied for the graph of Fig. 1 yields a rough estimate: 

~ (-^) 6 / dpf(p,M) [m soft ] 2 exp[-87r 2 //(p)], (20) 

where [m so f t ] 2 = W^rrix^B* and / is a dimensionless function of the instanton size p and 
the masses {M} of quantum fluctuations in the graph. If pM < 1 for all fluctuations, / 
would be of order unity, while it is suppressed by some powers of 1/pM when pM ^> 1 
for some of fluctuations. We already noted that the Higgs masses are of order Aqcd- For 
p <C A.qq D , the negative beta function of supersymmetric SU(3) theory with Nf = 6 is not 
negligible, implying exp[— 8tt 2 / g 2 (p)} leads to a significant suppression. For larger instantons 
with p ^> Aqq D , although exp[— 8n 2 /g 2 (p)} is roughly a constant since the beta function is 
small enough to approach to the fixed point |L3| , / is suppressed by the large Higgs masses. 



The above arguments imply that the dominant contribution is from instantons with size 
p ~ Aq CD . We thus have 

V a ~ e ia ^(-^)V* 3 det(A M A d )[m soft ] 2 A Q ^ D + h.c. (21) 

where [m so f t ] 2 is given in Eq. (18). 

In the above, we have computed the axion potential for the MSSM with Aqcd ^ m soft- I n 
fact, the MSSM can be considered as a rather special case since it contains a supersymmetric 
dimensionful parameter p. To see what happens in models without such a parameter, let us 
consider the next minimal supersymmetric standard model (NSSM) including an additional 
gauge singlet S with the superpotential 

W = \ u H u Qu c + X d H d Qd c + MSHuHa + X 2 S 3 . (22) 

As the MSSM, one can use the supersymmetry and also the selection rules of 

Gnssm = SU(3) Q x SU(3) u c x SU(3) d c x U(1) A x U(l) x x U(l) x > x U(1) R , 
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to constrain the effective Kahler potential K e s leading to the axion potential. (For the 
quantum numbers of Gnssm, see Table 2.) We then find K eS is proportional to 

e- Y det(Z u Z d )ZfD 2 Z 2 F( V ), (23) 

where the new spurion superfields Z\ and Z 2 are defined as 

= (1 + A0 2 )Ai, Z 2 = (1 + A6 2 )X 2 , 

and D 2 = D a D a denotes the supercovariant derivative. Although D 2 is applied to Z 2 in 
the above example, it can be applied to other spurions also. For a rough estimate, one may 



consider Fig. 2. Again applying the dimensional analysis rule given in Ref. |15| , we find 

1 



v 16vr 

where again [m so f t ] 2 is given in Eq. (18). 



Ksoft] Aqod + h.c, (24) 



III. RELAXATION OF THE AXION MISALIGNMENT 

In the previous section, we have estimated the axion potential by distinguishing the two 
cases: (i) Aqcd < "isoft and (ii) Aqcd 3> "isoft- Based on this, in this section we examine 
whether the axion misalignment can be relaxed to a small value by a stronger QCD in the 
early universe. To proceed, let us note that in generic supergravity models some low energy 
parameters other than Aqcd are determined also by the VEVs of some fields. In the early 
universe, such parameters may take values which differ from the present ones. However, in 
order for the early axion VEV to coincide with the present VEV, one should require that 
complex parameters which affect the axion VEV have almost the same values as the present 
ones. For the MSSM, 

(9 c s)mssm = 9qcd + arg(det(A u A d )) + 3arg(m 1/2 ) - 3arg(/xB), (25) 

while for the NSSM we have 
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(Oes)NSSM = Oqcd + arg(det(A u A d )) + 3arg(mi/ 2 ) - 2arg(A) - arg(A 1 A^). (26) 

This leads us to assume that in the early universe all Yukawa couplings, /i, and complex soft 
parameters (A, B, mi/ 2 ) have the same values as the present ones in order that the early 
axion VEV can coincide with the present VEV. 

Clearly the energy density in the early universe contains the dilaton potential energy 
V((f)). At present with <p = <fio, we have V(4>o) = 0. However, in the early universe with 
7^ 0o if the dilaton potential is not flat, which is believed to be the case, V(<p) would have 
a nonzero positive value. This means that raising up Aqcd to raise up the axion mass m a 
raises up also the energy density and thus the expansion rate H. In supersymmetric models, 
it is convenient to parameterize this dilaton potential energy in the early universe as follows: 

V(<P) = Cml {t M 2 m , (27) 

where C is a dimensionless coefficient and M m is the messenger scale of SUSY breaking 
above which the soft SUSY breaking is not operative anymore and thus a precise SUSY 
cancellation takes place. To be more concrete, here we specify m soft as the size of soft SUSY 
breaking for supermultiplets in the supersymmetric standard model sector. 

Since it is relevant for our later discussion, let us estimate the size of C. First of all, if the 
dilaton potential is generated directly by the SUSY breaking dynamics, V(4>) is expected 
to be of order |F| 2 where F denotes the auxiliary components of generic fields in SUSY 
breaking sector. For hidden sector models in which this SUSY breaking is transmitted by 



supergravity interactions, we have M m ~ Mpi and m so f t ~ \F\/Mp\ [fL7|| , and thus C is of 



order unity. For visible sector models in which the SUSY breaking is transmitted by gauge 



interactions [f[q| , the messenger scale corresponds to the scale of dynamical SUSY breaking 
and thus M m ~ l-^l 1 ^ 2 , while the soft breaking in the supersymmteric standard model sector 
is radiatively generated as m so f t ~ (^) n |-^| 1//2 where n is a model-dependent positive integer. 
This implies that C is of order (47r/a) 2ri for visible sector models when the dilaton potential 
energy directly induced by the SUSY breaking dynamics is parameterized as Eq. (27). In 
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summary, if the dilaton potential is generated directly by the SUSY breaking dynamics, C 
would be of order unity or bigger by (47r/a) 2n . 

Even in the case that the dilaton potential is not induced directly by the SUSY breaking 
dynamics, once SUSY is broken, it is generated in general by higher order loop effects. For 
instance, loops of colored particles would induce a dilaton-dependent, i.e. QCD coupling- 
dependent, contribution to the vacuum energy density. Such loop effects are quadratic in 
both the mass splitting in supermultiplets and the messenger scale M m which corresponds 
to the cutoff scale, and thus can be written as Eq. (27). (Throughout this paper, the 
dilaton is minimally defined as a field whose VEV determines the QCD coupling. Of course 
this dilaton can affect other gauge couplings. Our whole discussion will be valid also for 
such general case.) In the hidden sector models, the mass splitting is of order m so f t which 
is independent of the QCD coupling at the leading order. As a result, a QCD coupling- 
dependence would not appear at one loop order, but does appear at two loops fL9|. This 



implies that, due to radiative effects, C cannot be significantly smaller than (jg^-) 2 in the 
hidden sector models. In visible sector models, things are a bit more complicated. Typical 
visible sector models include a so-called messenger sector which contains a vector-like quark 
and lepton superfields [|TSJ. By the aid of messegner U(l) gauge interaction, the SUSY 
breaking dynamics gives rise to a mass splitting d~M in the messenger sector. Subsequent 
radiative effects of the standard model gauge interactions then lead to the mass splitting 
^soft — -f~8M in the supersymetric standard model sector. Since SM is independent of the 
QCD coupling at leading order, loops of the messenger sector particles will induce a QCD 
coupling- dependent vacuum energy density again at two loop order. The corresponding 
dilaton potential energy is of order (j^) 2 (5 M) 2 which is of order m^M^. Thus in 
visible sector models, again due to radiative effects, C cannot be significantly smaller than 
order unity. 

With the dilaton potential energy given in Eq. (27) in the early universe, the Hubble 
expansion rate is given by 
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H ~ C^m^MmMpy 1 , (28) 

where Mpi = 2.44 x 10 18 GeV. Let us now examine whether the two conditions of Eq. (2) 
for the relaxation of the axion misalignment can be satisfied. We consider first the case 
that Aqcd < fisoft without any light quark. One can easily confirm that, as long as the 
complex parameters contributing to # e ff are unchanged, the axion VEV determined by the 
early axion potential Eq. (4), i.e. 6 in , coincides with 8 e g. In the previous section, we have 
noted that all quark masses would become larger than Aqcd if both of the Higgs soft masses 
m\ and m 2 , of Eq. (6) receives a significantly large negative contribution from the energy 
density in the early universe. This would make the Higgs VEVs slide down toward the flat 
direction \H U \ = \Hd\ ^> m so { t . (Note that a large positive contribution leads to vanishing 
Higgs VEVs.) However, as was noted in Ref. | T2]] , such a negative contribution arises usually 



through supergravity interactions of the form -^r J d 2 9d 2 9 (p(p*XX* where x is an inflaton-like 



pi 

field yielding an energy density of order H 2 Mp Y . For the above operators significantly alter 
the Higgs soft mass, one needs H to be comparable to m so f t . Obviously then, we obtain 
m a -C H since Aqcd < ^soft <C f a in this case. 

Let us consider the next case that Aq C d is still less than m soft but now there exists a 
light quark. The axion potential in this case is given by Eq. (5). Using this axion potential 
and Eq. (28), we then find 

^- 5C \ . —jWjKW 2 -) U^J KW) ' (29) 

where all numbers in the brackets denote the energy scales in GeV unit. The messenger scale 
M m can be as low as 10 5 GeV in the visible sector models for SUSY breaking. However, 
in other type of models including the popular hidden sector models, M m is typically much 
bigger than 10 5 GeV. As was mentioned, complex parameters contributing to 9 e g in the early 
universe are required to have the values which are the same as the present one in order for 
the early axion VEV (^ n ) to coincide with the present VEV (0 e ff)- Real soft scalar masses 
m 2 escape from this requirement. They might be significantly bigger than the present values 
of order 10 2 GeV due to contributions from the radiation or other forms of energy density. 
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Note that if the early universe does not carry quantum numbers which are carried by some 
complex parameters, a nonzero energy density in such an early universe would enhance mf 
without affecting those complex parameters. However, if the enhanced soft mass-squared of 
Higgs fields are significantly bigger than the unaffected Bfi ~ 10 2 GeV, it makes the Higgs 
VEVs to vanish. This leads to m q = and thus dramatically suppresses the axion potential. 

Taking into account the points discussed above, it is easy to note that, for a given value 
of f a , the ratio m a /H becomes maximal for Aqcd comparable to m so ft which is of order 10 2 
GeV. Then m a would be large enough to push the axion field toward the minimum if f a is in 
the range: fa < (5C- 1 / 2 ) x (10 5 /M m ) x (4 x 10 12 ) GeV. Clearly in the hidden sector models 
with M m = M P1 , this range of f a does not overlap with the interesting range /„>4x 10 12 
GeV. In visible sector models, the messenger scale M m corresponds to the scale of dynamical 
SUSY breaking and thus may be as low as 10 5 GeV. However as we have discussed, C is at 
least of order unity in visible sector models and thus the axion mass is not large enough to 
relax the axion misalignment for / a >4x 10 12 GeV. 

So far, we have argued that raising Aq C d up to the order of m soft or below (while keeping 
the complex parameters that contribute to 9 e g unchanged) is not useful for relaxing the axion 
misalignment when f a 3> 4 x 10 12 GeV. One might expect that raising Aqcd further up to far 
above m so ft can leads to an axion mass m a > H. However, the axion potentials Eq. (21) and 
Eq. (24) computed in the previous section indicate that the axion mass is highly suppressed 
for Aqcd 3> "Wt- This is mainly because the Higgs VEVs and the condensates of the 
quarks, squarks, and the gluinos do not become of order Aq C d, but they all vanish when 
Aqcd 3> ""isoft- As a result, instanton amplitudes for the axion potential are suppressed by 
the powers of small Yukawa couplings and also of the small mass parameters ji and m so f t . 

Before discussing the size of the axion mass, let us briefly discuss the axion VEV. As can 
be noticed, the axion VEV determined by the axion potentials Eq. (21) and Eq. (24) differ 
from the 9 e s of Eq. (25) and Eq. (26) even when all complex parameters contributing 9 e g in 
the early universe have the same values as the present one. In such a case, the misalignment 
angle 59 = (a/f a ) - 9 cS for the MSSM is given by 
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{S6) MSS m = N ia ig(m 1/2 A*) + N 2 aig(m 1/2 B*), (30) 
where Aq and N 2 are appropriate integers of order unity. Similarly for the NSSM, we have 

(S9) NSSM = Narg(m 1/2 A*) (31) 
for an integer N of order unity. Although can be nonzero, the above misalgnment angle is 



severely constrained by the neutron electric dipole moment p0 



as 



59 < 1(T 2 ~ 1(T 3 . (32) 

This would be small enough to raise the cosmological bound on f a up to the grand unification 
scale or the Planck scale. However as we mentioned, the axion mass becomes too small to 
relax down the axion misalignment to a small value of order 89. 
Again using the axion potentials Eq. (21) and Eq. (24), we find 

( m '\ ~ in-"r-V ( i x 1( H f 10 " \( " Y' 2 ( 102 V' 2 

{ it )mssm " 10 c v^r) \m) iw-) fej • 



1/2 

) , (33) 



where again the numbers in the brakets denote energy scales in GeV unit. The above results 
obviously indicate that m a H . Although only two models are explicitly considered for the 
case Aqcd 3> ^soft, the huge suppression of the axion potential seems to be quite generic. In 
particular, adding more colored particles leads to a further suppression, and thus not helpful 
at all. We thus conclude that the axion misalignment cannot be relaxed down to a small 
value even when Aqcd 3> m so f t . In summary, the analysis made in this section indicates 
that a stronger QCD in the early universe is not useful for relaxing the axion misalignment. 



IV. CONCLUSION 



In this paper, we have examined whether the axion misalignment can be relaxed down 
to a small value by a stronger QCD in the early universe. This would allow the axion 
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scale to be of order the grand unification scale or the Planck scale without any cosmological 
difficulty. We discussed in somewhat detail the axion potentials in the early universe, in 
particular for the case that Aqcd ^soft- Taking into account the dilaton potential energy 
associated with a stronger QCD, our analysis indicates that the two conditions m a > H 
and (a/ f a ) = 6> cff (up to a small misalignment of order 1CT 3 ~ 1CT 2 ) cannot be satisfied 
simultaneously. We thus conclude that a stronger QCD in the early universe is not useful 
for raising the cosmological upper bound of the axion scale. 
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TABLES 



TABLE I. Quantum numbers of superflelds and spurions in MSSM 





U(1)a 


U(l)x 


U(1)r 
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i 





1 


u c , d c 
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-1 
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H u , H d 
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e~ Y 


12 


-6 


6 




-2 













-2 


2 


d 2 e 








-2 



TABLE II. Quantum numbers of superfields and spurions in NSSM 





U(1)a 


U(l)x 


U(l)x> 


U(1)r 
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1 


u c , d c 
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-1 





1 


H u , Hd 





1 








S 








1 





e~ Y 


12 


-6 





6 


Z u , H d 


-2 











Zi 





-2 


-1 
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z 2 








-3 


2 


d 2 e 











-2 
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FIGURES 

FIG. 1. Instanton graph for the axion potential Eq. (21) of the MSSM. The solid lines with 
and without waves around the instanton denote the gluino and quark modes, respectively, while 
the dotted lines are the Higgs and squarks fluctuations. The dark blobs represent the insertions 
of complex couplings which are explicitly written in the graph. The vertices not marked with 
couplings are the QCD gauge couplings. 

FIG. 2. Instanton graph for the axion potential Eq. (24) of the NSSM. Again the solid lines 
are for fermion modes , the dotted lines for boson fluctuations, and the dark blobs for the inserted 
complex couplings. 



23 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-ph/9606372vl 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-ph/9606372vl 



